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Abstract. The maximum integer skew-symmetric flow problem (MSFP) generalizes both the 
maximum flow and maximum matching problems. It was introduced by Tutte J8j in terms of self- 
conjugate flows in antisymmetrical digraphs. He showed that for these objects there are natural 
analogs of classical theoretical results on usual network flows, such as the flow decomposition, 
augmenting path, and max-flow min-cut theorems. We give unified and shorter proofs for those 
theoretical results. 

We then extend to MSFP the shortest augmenting path method of Edmonds and Karp [Jj 
and the blocking flow method of Dinits [Ij, obtaining algorithms with similar time bounds in 
general case. Moreover, in the cases of unit arc capacities and unit "node capacities" the block- 
ing skew-symmetric flow algorithm has time bounds similar to those established in |S1 for 
Dinits' algorithm. In particular, this implies an algorithm for finding a maximum matching in a 
nonbipartite graph in 0{y/nm) time, which matches the time bound for the algorithm of Micali 
and Vazirani 25 . Finally, extending a clique compression technique of Feder and Motwani [2] 
to particular skew-symmetric graphs, we speed up the implied maximum matching algorithm to 
run in 0(y^mlog(n^/m)/logn) time, improving the best known bound for dense nonbipartite 
graphs. 

Also other theoretical and algorithmic results on skew-symmetric flows and their applications 
are presented. 
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1 Introduction 



By a skew- symmetric flow we mean a flow in a skew-symmetric directed graph which takes equal 
values on any pair of "skew-symmetric" arcs. This is a synonym of Tutte's self- conjugate flow in 
an antisymmetrical digraph [2Hj. This paper is devoted to the maximum integer skew-symmetric 
flow problem, or, briefly, the maximum IS-flow problem. We study combinatorial properties of 
this problem and develop fast algorithms for it. 

A well-known fact jlUj is that the bipartite matching problem can be viewed as a special case 
of the maximum flow problem. The combinatorial structure of nonbipartite matchings revealed 
by Edmonds ^ involves blossoms and is more complicated than the structure of flows. This 
phenomenon explains, to some extent, why general matching algorithms are typically more in- 
tricate relative to flow algorithms. The maximum IS-flow problem is a generalization of both 
the maximum flow and maximum matching (or b-matching) problems. Moreover, this general- 
ization appears to be well-grounded for two reasons. First, the basic combinatorial and linear 
programming theorems for usual flows have natural counterparts for IS-flows. Second, when solv- 
ing problems on IS-flows, one can use intuition, ideas and technical tools well-understood for usual 
flows, so that the implied algorithms for matchings become more comprehensible. 

As the maximum flow problem is related to certain path problems, the maximum IS-flow 
problem is related to certain problems on so-called regular paths in skew-symmetric graphs. We 
use some theoretical and algorithmic results on the regular reachability and shortest regular path 
problems from |16j . 

Tutte j28j originated a mini-theory of IS-flows (in our terms) to bridge theoretical results 
on matchings and their generalizations (b- factors, b-matchings, degree-constrained subgraphs, 
Gallai's track packings, and etc.) and results on usual flows. This theory parallels Ford and 
Fulkerson's flow theory jlOj and includes as basic results the decomposition, augmenting path, 
and max- flow min-cut theorems. Subsequently, some of those results were re-examined in different, 
but equivalent, terms by other authors, e.g., in [HI IT^ 1^ . 

Recall that the flow decomposition theorem says that a flow can be decomposed into a col- 
lection of source-to-sink paths and cycles. The augmenting path theorem says that a flow is 
maximum if and only if it admits no augmenting path. The max- flow min-cut theorem says that 
the maximum flow value is equal to the minimum cut capacity. Their skew-symmetric analogs 
are, respectively, that an IS-flow can be decomposed into a collection of pairs of symmetric source- 
to-sink paths and pairs of symmetric cycles, that an IS-flow is maximum if and only if it admits 
no regular augmenting path, and that the maximum IS-flow value is equal to the minimum odd- 
barrier capacity. We give unified and shorter proofs for these skew-symmetric flow theorems. 

There is a relationship between skew-symmetric flows and bidirected flows introduced by Ed- 
monds and Johnson [Hj in their combinatorial study of a natural class of integer linear programs 
generalizing usual flow and matching problems. In particular, they established a linear program- 
ming description for integer bidirected flows. We finish the theoretical part by showing how to 
obtain a linear programming description for maximum IS-flows directly, using the max-IS-flow 
min-barrier theorem. 
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The second, larger, part of this paper is devoted to efficient methods to solve the maximum IS- 
flow problem (briefly, MSFP) in general and special cases, based on the theoretical ground given 
in the first part. First of all we explain how to adapt the idea of Anstee's elegant methods 012] 
for b-matchings in which standard flow algorithms are used to construct an optimal half-integer 
solution and then, after rounding, the "good pre-solution" is transformed into an optimal b- 
matching by solving 0(n) certain path problems. We devise an 0(M(n, m) +nm)-time algorithm 
for MSFP in a similar fashion, using a regular reachability algorithm with linear complexity to 
improve a good pre-solution. Hereinafter n and m denote the numbers of nodes and arcs of the 
input graph, respectively, and M(n, m) is the time needed to find an integer maximum flow in a 
usual network with n nodes and m arcs. Without loss of generality, we assume n = 0{m). 

The next approach is the core of this paper. The purpose is to extend to MSFP the well-known 
shortest augmenting path algorithm of Edmonds and Karp [7j with complexity 0{nm?), and its 
improved version, the blocking flow algorithm of Dinits [1] with complexity 0{n^m), so as to 
preserve the complexity bounds. Recall that the blocking flow algorithm consists of 0{n) phases, 
each flnding a blocking flow in a layered network representing the union of currently shortest 
augmenting paths. We introduce concepts of shortest blocking and totally blocking IS-flows and 
show that an optimal solution to MSFP is also constructed in 0{n) phases, each finding a shortest 
blocking IS-fiow in the residual skew-symmetric network. In its turn a phase is reduced to finding 
a totally blocking IS-fiow in an acyclic (though not necessarily layered) skew-symmetric network. 

The crucial point is to perform the latter task in time comparable with the phase time in 
Dinits' algorithm (which is 0{nm) in general case). We reduce it to a certain auxiliary problem 
in a usual acyclic digraph. A fast algorithm for this problem provides the desired time bound for 
a phase. 

The complexity of our blocking IS-fiow algorithm remains comparable with that of Dinits' 
algorithm in important special cases where both the number of phases and the phase time signif- 
icantly decrease. More precisely, Dinits' algorithm applied to the maximum matching problem in 
a bipartite graph runs in 0{y/nm) time jl8l l2Uj . Extending that result, it was shown in [SI |^ 
that for arbitrary nonnegative integer capacities, Dinits' algorithm has 0(min{n, \/A}) phases 
and each phase runs in O (min{ nm, m + A}) time, where A is the sum of transit capacities of 
inner nodes. Here the transit capacity of a node (briefly, the node capacity) is the maximum flow 
value that can be pushed through this node. We show that both bounds remain valid for the 
blocking IS-fiow algorithm. 

When the network has unit arc capacities (resp. unit inner node capacities), the number of 
phases turns into 0{^/m) (resp. 0{^/n))■, in both cases the phase time turns into 0{m). The 
crucial auxiliary problem (that we are able to solve in linear time for unit arc capacities) becomes 
the following maximal balanced path-set problem: 

MBP: Given an acyclic digraph in which one sink and an even set of sources partitioned into pairs 
are distinguished, find an (inclusion-wise) maximal set of pairwise arc-disjoint paths from 
sources to the sink such that for each pair {z,z'} of sources in the partition, the number of 
paths from z is equal to that from z' . 
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As a consequence, the implied algorithm solves the maximum matching problem in a general 
graph in the same time, 0{^/nm), as the algorithm of Micali and Vazirani [251 129j (cf. [HIE]) and 
solves the b-factor or maximum degree-constrained subgraph problem in 0(m^/^) time, similarly 
to Gabow . The logical structure of our algorithm differs from that of and sophisticated 
data structures (incremental trees for set union ^5]) are used only in the regular reachability and 
shortest regular path algorithms of linear complexity from 16^ (applied as black-box subroutines) 
and once in the algorithm for MBP. 

Finally, we show that a clique compression technique of Feder and Motwani P| can be extended 
to certain skew-symmetric graphs. As a result, our maximum matching algorithm in a general 
graph is speeded up to run in 0(\/rem log (n^/m)/ log n) time. This matches the best bound for 
bipartite matching [S]. 

Fremuth-Paeger and Jungnickel developed an algorithm for MSFP (stated in terms of 
"balanced flows") which combines Dinits' approach with ideas and tools from j25[ I29j : it runs in 
0{nm?) time for general capacities and in time slightly slower than 0{^/nm) in the nonbipartite 
matching case. 

This paper is organized as follows. Basic definitions and facts are given in Section [21 Sec- 
tions 01 and 01 contain theoretical results on combinatorial and linear programming aspects of 
IS-flows, respectively. Section describes Anstee's type algorithm for MSFP. The shortest reg- 
ular augmenting path algorithm and a high level description of the blocking IS-flow method are 
developed in Sectional Section [3 gives a short review on facts and algorithms in jl6j for regular 
path problems. Using it, Section |H1 explains the idea of implementation of a phase in the block- 
ing IS-flow method. It also bounds the number of phases for special skew-symmetric networks. 
Section El completes the description of the blocking IS-fiow algorithm by reducing the problem of 
finding a totally blocking IS-flow in an acyclic skew-symmetric network to the above-mentioned 
auxiliary problem in a usual acyclic digraph and devising a fast procedure to solve the later. The 
concluding Sectional discusses implications for matchings and their generalizations, and explains 
how to speed up the implied maximum matching algorithm by use of the clique compression. 

This paper is self-contained up to several quotations from ^Sl- Main results presented in this 
paper were announced in extended abstract jTH]. Subsequently, the authors found a flaw in the 
original fast implementation of a phase in the blocking IS-flow method. It was corrected in a 
revised version of this paper (circulated in 2001) where problem MBP and its weighted analog 
were introduced and efficiently solved, whereas the original version (identical to preprint |17j ) 
embraced only the content of Sections (HHSl 

2 Preliminaries 

By a skew- symmetric graph we mean a digraph G = (V, E) with a mapping (involution) a of 

V \J E onto itself such that: (i) for each x G V U E, cr(x) 7^ x and a{(T{x)) = x; (ii) for each 

V £ V, a{v) G V; and (iii) for each a = {v,w) G E, a{a) = {a{w),a{v)). Although parallel arcs 
are allowed in G, an arc leaving a node x and entering a node y is denoted by (x, y) when it is 
not confusing. We assume that a is fixed (when there are several such mappings) and explicitly 
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included in the description of G. The node (arc) (t{x) is called symmetric to a node (arc) x (using, 
for brevity, the term symmetric rather than skew-symmetric). Symmetric objects are also called 
mates, and we usually use notation with primes for mates: x' denotes the mate (j{x) of an element 
X. Note that G can contain an arc a from a node v to its mate v'\ then a' is also an arc from v 
to v' . 

Unless mentioned otherwise, when talking about paths (cycles), we mean directed paths (cy- 
cles). The symmetry a is extended in a natural way to paths, subgraphs, and other objects in G; 
e.g., two paths (cycles) are symmetric if the elements of one of them are symmetric to those of 
the other and go in the reverse order. Note that G cannot contain self-symmetric paths or cycles. 
Indeed, if P = (xq, ai, xi, . . . , a^, x^) is such a path (cycle), choose arcs Oj and aj such that i < j, 
Uj = cr(aj) and j — i is minimum. Then j > i + 1 (as j = i would imply cr{ai) = ai and j = i + 1 
would imply cr(xj) = Xj-i = Xi). Now cr(aj+i) = Oj-i contradicts the minimality of j — i. 

We call a function h on E symmetric if h{a) = h(a') for all a £ E. 

A skew- symmetric network is a quadruple = (G, a, u, s) consisting of a skew-symmetric 
graph G = {V, E) with symmetry a, a nonnegative integer-valued symmetric function u (of arc 
capacities) on E, and a source s £ V. The mate s' of s is the sink of N. A flow in A*" is a function 
f : E ^ M-l- satisfying the capacity constraints 

/(a) < u{a) for all a £ E 

and the conservation constraints 

div/(x) := ^ f{x,y)- ^ /(y,x) = ior all x £ V - {s, s'}. 

{x,y)eE iy,x)eE 

The value divj-(s) is called the value of / and denoted by |/|; we usually assume that |/| > 0. 
Now IS-flow abbreviates integer symmetric flow, the main object that we study in this paper. 
The maximum IS-flow problem (MSFP) is to find an IS-flow of maximum value in N . 

The integrality requirement is important: if we do not require / to be integral, then for any 
integer flow / in N, the flow /', defined by /'(a) := (/(a) + f{a'))/2 for a G is a flow of 
the same value as /, which is symmetric but not necessarily integral. Therefore, the fractional 
skew-symmetric flow problem is equivalent to the ordinary flow problem. 

Note that, given a digraph D = {V{D),A(D)) with two specified nodes p and q and non- 
negative integer capacities of the arcs, we can construct a skew-symmetric graph G by taking a 
disjoint copy D' of D with all arcs reversed, adding two extra nodes s and s' , and adding four arcs 
(s,p), {s, q'), (g, s'), (p', s') of infinite capacity, where p', q' are the copies of p, q in D', respectively. 
Then there is a natural one-to-one correspondence between integer flows from p to q m D and 
the IS-flows from s to s' in G. This shows that MSFP generalizes the classical (integer) max-flow 
problem. 

Remark. Sometimes it is useful to consider a sharper version of MSFP in which double-sided 
capacity constraints £{a) < f{a) < u{a), a £ E, are imposed, where i,u : E ^ Z+ and £ < u {prob- 
lem DMSFP). Similarly to the max-flow problem with upper and lower capacities |1U| . DMSFP 
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is reduced to MSFP in the skew-symmetric network A'^' obtained from N by subdividing each 
arc a = {x,y) into three arcs {x,v), {v,w), {w,y) with (upper) capacities u{a),u{a) — £{a),u{a), 
respectively, and adding extra arcs {s,w) and {v,s') with capacity i{a) each. It is not difficult 
to show (e.g., using Theorem 13. 5|) that DMSFP has a solution if and only if all extra arcs are 
saturated by a maximum IS-flow /' for N' , and in this case /' induces a maximum IS-flow for 
in a natural way. For details, see jTTj . 

In our study of IS-flows we rely on results for regular paths in skew-symmetric graphs. A 
regular path, or an r-path, is a path in G that does not contain a pair of symmetric arcs. Similarly, 
an r-cycle is a cycle that does not contain a pair of symmetric arcs. The r-reachability problem 
(RP) is to find an r-path from s to s' or a proof that there is none. Given a symmetric function 
of arc lengths, the shortest r-path problem (SRP) is to find a minimum length r-path from s to s' 
or a proof that there is none. 

A criterion for the existence of a regular s to s' path is less trivial than that for the usual path 
reachability; it involves so-called barriers. We say that 

B = {A-Xi,...,Xk) 

is an s -barrier if the following conditions hold. 

(Bl) A,Xi, . . . , Xk are pairwise disjoint subsets of V, and s £ A. 
(B2) For A' = a{A), AnA' = i!}. 

(B3) For i = 1, . . . ,k, Xi IS self-symmetric, i.e., o"(Xj) = Xi. 
(B4) For i = 1, . . . ,k, there is a unique arc, e*, from A to Xj. 
(B5) For i,j = l,...,k and i ^ j, no arc connects Xi and Xj. 

(B6) For M := y - (A U A' U Xi U . . . U Xj,) and i = 1, . . . , fc, no arc connects Xi and M. 
(B7) No arc goes from Ato A' U M. 

(Note that arcs from A' to A, from Xi to A, and from M to ^ are possible.) Figure ^ illustrates 
the definition. Tutte proved the following (see also pilTB]). 

Theorem 2.1 j2H] There is an r-path from s to s' if and only if there is no s-barrier. 

This criterion will be used in Section to obtain an analog of the max-flow min-cut theorem 
for IS-flows. RP is efficiently solvable. 

Theorem 2.2 [51 116j The r-reachability problem in G can be solved in 0{m) time. 

The methods for the maximum IS-flow problem that we develop apply, as a subroutine, the 
r-reachability algorithm of linear complexity from |16j . which flnds either a regular s to s' path 
or an s-barrier. Another ingredient used in our methods is the shortest r-path algorithm for the 
case of nonnegative symmetric lengths, which runs in 0{m logn) time in general, and in 0{m) 
time for all- unit lengths ^SJ- The necessary results on RP and SRP are outlined in Section [7J 
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Figure 1: A barrier 

In the rest of this paper, a and s wiU denote the symmetry map and the source, respectively, 
regardless of the network in question, which will allow us to use the shorter notation (G, u) for a 
network (G,a,u,s). Given a simple path P, the number of arcs on P is denoted by \P\ and the 
incidence vector of its arc set in is denoted by i i-e-; X^(fl) = 1 if a is an arc of P, and 
otherwise. 

2.1 Relationships to Matchings 

Given an undirected graph G' = {V , E'), a matching is a subset M Q E' such that no two edges 
of M have a common endnode. The maximum matching problem is to find a matching M whose 
cardinality |M| is as large as possible. 

There are well-known generalizations of matchings; for a survey see j231 1241 OS] . Let uo,u : 
E' U {oo} and bo,b : V be functions such that b^ < b and uq < u. A {uq^u)- 

capacitated (bo, b) -matching is a function h : E' ^ Z+ satisfying the capacity constraints 

uo{e) < h{e) < u{e) for all e G E\ 

and the supply constraints 

bo{v) < h{e) < b{v) for all v G V . 

e={v,w}£E' 

The value of h is defined to be h{E'). Hereinafter, for a numerical function g' on a set S and a 
subset S' C S, g{S') denotes J2eeS' di^)- Popular special cases are: a u- capacitated b-matching 
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Matching Skew-symmetric flow 

Figure 2: Reduction example for maximum cardinality matching 

(when 60 = 0); a degree- constrained subgraph (when n = 1); a perfect b-matching (when u = 00 
and bo = b); a b- factor (when u = 1 and bo = b). In these cases one assigns uq = 0. Typically, in 
unweighted versions, one is asked for maximizing the value of h (in the former two cases) or for 
finding a feasible h (in the latter two cases). 

The general maximum (uq, 'u)-capacitated (60, 6)-matching problem is reduced to the maximum 
IS-flow problem (MSFP or DMSFP, depending on whether both uq, bo are zero functions or not) 
without increasing the problem size by more than a constant factor. The construction of the 
corresponding capacitated skew-symmetric graph G = {V, E) is straightforward (and close to 
that in |2H]): 

(i) for each v £ V, V contains two symmetric nodes vi and V2', 

(ii) also V contains two additional symmetric nodes s and s' (the source and the sink); 

(iii) for each e = {v,w} € E', E contains two symmetric arcs {vi,W2) and (1^1,^2) with lower 
capacity uo{e) and upper capacity u{e); 

(iv) for each v G V, E contains two symmetric arcs {s,vi) and {v2,s') with lower capacity 
bo{v) and upper capacity b{v). 

There is a natural one-to-one correspondence between the (uq, ii)-capacitated {bo, 6)-matchings 
h in G' and the IS-flows / from s to s' in G, and the value of / is twice the value of h. Figure [21 
illustrates the correspondence for matchings. 

In case of the b-factor or degree-constrained subgraph problem, one may assume that b does 
not exceed the node degree function of G. Therefore, one can make a further reduction to MSFP 
in a network with Od^''!) nodes, Od-E'l) arcs, and unit arc capacities (by getting rid of lower 
capacities as in the Remark above and then splitting each arc a with capacity q{a) > 1 into q{a) 
parallel arcs of capacity one). In Section we compare the time bounds of our methods for 
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MSFP applied to the matching problem and its generalizations above with known bounds for 
these problems. 

Edmonds and Johnson ^ studied the class of integer linear programs in which the constraint 
matrix entries are integers between -2 and +2 and the sum of absolute values of entries in each 
column (without including entries from the box constraints) does not exceed two. Such a problem 
is often stated in terms of bidirected graphs (for a survey, see [261 Chapter 36]). Recall that a 
bidirected graph H = {X, B) may contain, besides usual directed edges going from one node to 
another, edges directed from both of its endnodes, and to both of them. A particular problem on 
such an object is the maximum bidirected flow problem: given a capacity function c : B ^ Z_|_ and 
a terminal p £ X, find a function (biflow) g : B ^ Z+ maximizing the value divg{p). (Reasonable 
versions with more terminals are reduced to this one.) Here g < c and divg{x) = for all 
X £ X — {p}, where divg{x) is the total biflow on the edges directed from x minus the total biflow 
on the edges directed to x (a loop e at x contributes 0, 2g{e) or —2g{e)). 

The maximum IS-flow problem admits a linear time and space reduction to the maximum 
biflow problem (in fact, both are equivalent). More precisely, given an instance N = (G = 
{V, E), a, u, s) of MSFP, take a partition (X, X') of V such that X' = o'{X) and s £ X. For each 
pair {a, a'} of symmetric arcs in E and nodes x,y £ X, assign an edge from x to y if a or a' goes 
from X to y; an edge from both x,y if a or a' goes from x to C7(y); an edge to both x,y if a or a' 
goes from a{x) to y. This produces a bidirected graph H = {X, B). We set p := s and assign the 
capacity c(e) of each edge e G -B to be the capacity of the arc from which e is created. There is a 
one-to-one correspondence between the IS-flows in N and the billows in (H,c,p), and the values 
of corresponding flows are equal. A reverse reduction is also obvious. Using these reductions, one 
can try to derive results for IS-flows from corresponding results on biflows, and vice versa. In this 
paper we give direct proofs and algorithms for IS-flows. 

3 Mini-Theory of Skew- Symmetric Flows 

This section extends the classical flow decomposition, augmenting path, and max-flow min-cut 
theorems of Ford and Fulkerson ^U] to the skew-symmetric case. The support {e £ S : /(e) ^ 0} 
of a function / : S" ^ M is denoted by supp(/). 

Let h he a nonnegative integer symmetric function on the arcs of a skew-symmetric graph 
G = {V, E). A path (cycle) P in G is called h-regular if h{a) > for all arcs a of P and each arc 
a £ P such that a' £ P satisfies h{a) > 2. Clearly when h is all-unit on E, the sets of regular 
and /i-regular paths (cycles) are the same. We call an arc a of P ordinary if a' P and define 
the h- capacity 5h{P) of P to be the minimum of all values h{a) for ordinary arcs a on P and all 
values \h{a)/2\ for nonordinary arcs a on P. 

To state the symmetric flow decomposition theorem, consider an IS-flow / in a skew-symmetric 
network N = {G = (y, P),u). An IS-flow g m. N \s called elementary if it is representable as 
g = 5x^ + ^X^' 1 where P is a simple cycle or a simple path from s to s' or a simple path from 
s' to s, P' = cr{P), and (5 is a positive integer. Since g is feasible, P is u-regular and 5 < 5u{P)- 
We denote g by {P,P',5). By a symmetric decomposition of / we mean a set D of elementary 
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flows such that / = '^{g : g S D). The following symmetric decomposition theorem (see IISI ) 
slightly generalizes a result by Tutte that there exists a symmetric set of |/| paths from s to 
s' such that any arc a is contained in at most /(a) paths. 

Theorem 3.1 For an IS-flow f in G, there exists a symmetric decomposition consisting of at 
most m elementary flows. 

Proof. We build up an /-regular path F in G until this path contains a simple cycle P or a simple 
path P connecting s and s' . This will determine a member of the desired flow decomposition. 
Then we accordingly decrease / and repeat the process for the resulting IS-flow /', and so on 
until we obtain the zero flow. 

We start with F formed by a single arc a € supp(/). First we grow F forward. Let b = {v,w) 
be the last arc on the current (simple) path F. Suppose that w ^ s,s'. By the conservation for / 
at w, supp(/) must contain an arc q = {w, z). If q' is not on F or /(g) > 2, we add q to F. 

Suppose q' is on F and f[q) = 1. Let Fi be the part of F between w' and w. Then Fi contains 
at least one arc since w ^ w' . Suppose there is an arc q £ supp(/) leaving w and different from q. 
Then we can add g to F instead of q, forming a longer /-regular path. (Note that since the path 
F is simple, q' is not on F). Now suppose that such a g does not exist. Then exactly one unit of 
the flow / leaves w. Hence, exactly one unit of the flow / enters w, implying that b = (v, w) is 
the only arc entering w in supp(/), and that f{b) = 1. But a{d) also enters w, where d is the first 
arc on Fi. The fact that cr{d) ^ b (since Fi is /-regular) leads to a contradiction. 

Let {w, z) be the arc added to F. If z is not on F, then F is a simple /-regular path, and we 
continue growing F. If z is on F, we discover a simple /-regular cycle P. 

If F reaches s' or s, we start growing F backward from the initial arc a in a way similar to 
growing it forward. We stop when an /-regular cycle P is found or one of s, s' is reached. In the 
latter case P = F is either an /-regular path from s to s' or from s' to s, or an /-regular cycle 
(containing s or s'). 

Form the elementary flow g = {P, P' , 6) with S = Sf{P) and reduce / to /' := / — 5x^ — ^X^' ■ 
Since P is /-regular, 5 > 0. Moreover, there is a pair e, e' of symmetric arcs of P such that either 
/'(e) = f'{e') = or /'(e) = /'(e') = 1; we associate such a pair with g. In the former case e,e' 
vanish in the support of the new IS-flow /', while in the latter case e, e' can be used in further 
iterations of the decomposition process at most once. Therefore, each pair of arc mates of G is 
associated with at most two members of the constructed decomposition D, yielding \D\ < m. M 

The above proof gives a polynomial time algorithm for symmetric decomposition. Moreover, 
the above decomposition process can be easily implemented in 0{nm) time, which matches the 
complexity of standard decomposition algorithms for usual flows. 

The decomposition theorem and the fact that the network has no self-symmetric cycles imply 
the following useful property noticed by Tutte as well. 

Corollary 3.2 j2Hl For any self- symmetric set S V and any IS-flow in G, the total flow on 
the arcs entering S, as well as the total flow on the arcs leaving S, is even. 
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Remark. Another consequence of Theorem 13.11 is that one may assume that G has no arc 
entering s. Indeed, consider a maximum IS-flow / in G and a symmetric decomposition D of 
/. Putting together the elementary flows from s to s' in D, we obtain an IS-flow /' in G with 
l/'l ^ I/I) so /' is a maximum flow. Since /' uses no arc entering s or leaving s', deletion of all 
such arcs from G produces an equivalent problem in a skew-symmetric graph. 

Next we state a skew-symmetric version of the augmenting path theorem. It is convenient to 
consider the graph G^ = {V, E'^) formed by adding a reverse arc (y, x) to each arc {x, y) of G. For 
a G E'^, denotes the corresponding reverse arc. The symmetry a is extended to E'^ in a natural 
way. Given a (nonnegative integer) symmetric capacity function u on E and an IS-flow / on G, 
define the residual capacity Uf{a) of an arc a € E~^ to be u{a) — f{a) if a G i?, and f{a^) otherwise. 
An arc a G E~^ is called residual if Uf{a) > 0, and saturated otherwise. Given an IS-flow g in the 
network {G~^,Uf), we define the function f (B g on E by setting {f (B g){a) := f{a) + g{a) — g{a^). 
Clearly f (B g is a feasible IS-flow in {G, u) whose value is |/| -|- \g\. 

By an r- augmenting path for / we mean a n/-regular path from s to s' in G^ . If P is an 
r-augmenting path and if (5 G N does not exceed the uz-capacity of P, then we can push 5 units of 
flow through a (not necessarily directed) path in G corresponding to P and then 5 units through 
the path corresponding to P' . Formally, / is transformed into f (B g, where g is the elementary 
flow {P,P',6) in {G~^,Uf). Such an augmentation increases the value of / by 25. 

Theorem 3.3 |28j An IS-flow f is maximum if and only if there is no r-augmenting path. 

Proof. The direction that the existence of an r-augmenting path implies that / is not maximum 
is obvious in light of the above discussion. 

To see the other direction, suppose that / is not maximum, and let /* be a maximum IS- 
flow in G. For a e E define g{a) := f*{a) - f{a) and g{a^) := if f*{a) > /(a), while 
g{a^) := f{a) — f*{a) and g{a) := if /*(a) < /(a). One can see that g is a feasible symmetric 
flow in {G^,Uf). Take a symmetric decomposition D of g. Since \g\ = \ f*\ — \ f\ > 0, D has a 
member {P,P',6), where P is a uj-regular path from s to s'. Then P is an r-augmenting path 
for /. ■ 

In what follows we will use a simple construction which enables us to reduce the task of finding 
an r-augmenting path to the r-reachability problem. For a skew-symmetric network (H, h) , split 
each arc a = (x, y) of H into two parallel arcs ai and 02 from x to y (the first and second 
split-arcs generated by a). These arcs are endowed with the capacities [h]{ai) := \h{a)/2] and 
[/i](a2) := [h{a)/2\. Then delete all arcs with zero capacity [h]. The resulting capacitated graph 
is called the split-graph for {H,h) and denoted by S{H,h). The symmetry a is extended to the 
arcs of S{H, h) in a natural way, by defining a(ai) := {a{a))i for i = 1,2. 

For a path P in S{H, h), its image in H is denoted by io{P) (i.e., io{P) is obtained by replacing 
each arc Oj of P by the original arc a =: uj{ai)). It is easy to see that if P is regular, then uj{P) 
is /i-regular. Conversely, for any /i-regular path Q in H, there is a (possibly not unique) r-path 
P in S{H, h) such that lo{P) = Q. Indeed, replace each ordinary arc a of Q by the first split-arc 
ai (existing as h{a) > 1) and replace each pair a, a' of arc mates in Q by at, a'j for {i,j} = {1, 2} 
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(taking into account that h{a) = h{a') > 2). This gives the required r-path P. Thus, Theorem 13 .31 
admits the following reformulation in terms of split-graphs. 

Corollary 3.4 An IS-flow f in (G, u) is maximum if and only if there is no regular path from s 
to s' in S{G^ ,Uf). 

Finally, the classic max-flow min-cut theorem states that the maximum flow value is equal to 
the minimum cut capacity. A skew-symmetric version of this theorem involves a more complicated 
object which is close to an s-barrier occurring in the solvability criterion for the r-reachability 
problem given in Theorem 12.11 We say that B = {A; Xi, . . . , X^) is an odd s-barrier for (G, u) if 
the following conditions hold. 

(01) A,Xi,..., Xfc are pairwise disjoint subsets of V, and s £ A. 

(02) For A' = a{A), AnA' = i!}. 

(03) For i = 1, . . . ,k, Xi is self-symmetric, i.e., cr{Xi) = Xi. 

(04) For i = 1, . . . ,k, the total capacity u{A, Xi) of the arcs from A to Xi is odd. 

(05) For i,j = l,...,k and i ^ j, no positive capacity arc connects Xi and Xj. 

(06) For M := V — {A U A' L) Xi U . . . U X^) and i = 1, . . . , k, no positive capacity arc connects 
Xi and M. 

Compare with (B1)-(B7) in Sectional Define the capacity u{B) of B to be u{A, V — A) — k. Since 
the source is denoted by s throughout, we refer to an odd s-barrier as odd barrier. 
The following is the maximum IS-flow minimum barrier theorem. 

Theorem 3.5 |28| The maximum IS-flow value is equal to the minimum odd barrier capacity. 

Proof. To see that the capacity of an odd barrier B = (A; Xi, . . . , Xf^) is an upper bound on the 
value of an IS-flow /, consider a symmetric decomposition D of /. For each member g = (P, P' , 6) 
of D, where P is a path from s to s' , take the last arc a = {x,y) of the first path P such that 
X G A. If y £ A' , then the symmetric arc a' (which is in P') also goes from A to A' (by (02)), and 
therefore, g uses at least 25 units of the capacity of arcs from A to A' . Associate g with the pair 
a, a' . Now let y A' . Since y ^ A, y is either in y := M or in Y := Xi for some i. The choice 
of a and (01), (05), (06) imply that P leaves Y by an arc b from Y to A! . Then the symmetric 
arc b' (which is in P') goes from A to Y (since Y is self-symmetric), and therefore, g uses at 
least 2(5 units of the capacity ti(A, Y\ Associate g with the pair a, 6' (possibly a = b'). Note that 
at least one unit of each capacity u{A,Xi) is not used under the canonical way we associate the 
elementary s to s' flows of D with arcs from A to V — A (since u{A,Xi) is odd, by (04)). By 
these reasonings, |/| < u{B). 

Next we show that the two values in the theorem are equal. Let / be a maximum IS-flow. 
By Corollarv 13.41 the split-graph S = S{G^, uj) contains no s to s' r-path, so it must contain an 
s-barrier B = {A; Xi, . . . , Xk), by Theorem 12. 11 

Let e* be the (unique) arc from A to Xi in S (see (B4) in Section 2). By the construction of 
S, it follows that the residual capacity Uf of every arc from A to Xi in G^ is zero except for the 
arc a;(e*), whose residual capacity is one. Hence, 
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(i) if was formed by splitting an arc a £ E, then a goes from A to Xi, and /(a) = u{a) — 1; 

(ii) if e* was formed by splitting for a € E, then a goes from Xi to yl, and f{a) = 1; 

(iii) all arcs from A to Xi in G, except a in case (i), are saturated by /; 

(iv) all arcs from Xi to A in G, except a in case (ii), are free of flow. 

Furthermore, comparing arcs in S and G, we observe that: 

(v) property (B7) implies that the arcs from A to A' DM are saturated and the arcs from A' DM 
to A are free of flow; 

(vi) property (B5) implies (05) and (B6) implies (06). 

Properties (i)-(iv),(05),(06) together with Corollary 13.21 provide (04). So B is an odd s- 
barrier in G. We have |/| = f{A, V - A) - f{V - A,A) = u{A, V - A) - k {in view of (i)-(v)). 
Hence, |/| = u{B). ■ 

4 Integer and Linear Programming Formulations 

Although methods of solving MSFP developed in subsequent sections will not use explicitly linear 
programming aspects exhibited in this section, such aspects help to understand more about the 
structure of IS-flows. 

MSFP is stated as an integer program in a straightforward way. We use function rather than 
vector notation. For functions g,h on a set S, g ■ h denotes the inner product "^^^g g{x)h{x). 
Assuming that no arc of G enters the source s (see the Remark in the previous section), MSFP 
can be written as follows: 

maximize 1/1 = ^^^ ^ ^fis,v) subject to (1) 

/(a) > Va G ^ (2) 
/(a) < u{a) yaeE (3) 
-J2, + ,,p/("'^) = yv€V-{s,s'} (4) 

fia) - /(c7(a)) = yaGE (5) 
/(a) integer £ E (6) 

A linear programming formulation for MSFP is obtained by replacing the integrality condition 
(jHl) by linear constraints related to certain objects that we call odd fragments in G. The correctness 
of the resulting linear program will be shown by use of the max-min relation between IS-flow and 
odd barriers in Theorem 13.51 Alternatively, one can try to derive it from a linear programming 
characterization of integer bidirected flows in ^ (using the reduction as in Section [2J . 

An odd fragment is a pair p = {Vp, Up), where Vp is a self- symmetric set of nodes with s Vp, 
and Up is a subset of arcs entering Vp such that the total capacity u{Up) is odd. The characteristic 
function Xp of p is the function on E defined by 
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1 if a e UpUaiUp), 
Xpia) :={ -1 if a G 6{Vp) - {Up U a{Up)), (7) 
otherwise. 

Here d(Vp) is the set of arcs with one end in Vp and the other in V — Vp. We denote the set of 
odd fragments by 17. 

Let / be a (feasible) IS-flow and p ^ VL. By (O and the symmetry of u, we have f • xp ^ 
u{Up) + u{a{Up)) = 2u{Up). Moreover, / • Xp is at most 2u{Up) — 2; this immediately follows from 
Corollarv \^.2\ and the fact that u{Up) is odd. This gives new linear constraints for MSFP: 

f ■Xp< '^u{Up) - 2 for each p G 0. (8) 

Addition of these constraints enables us to drop off the symmetry constraints © and the inte- 
grality constraints © without changing the optimum value of the linear program. This fact is 
implied by the following theorem. 

Theorem 4.1 Every maximum IS-flow is an optimal solution to the linear program Q-©, ©• 

Proof. Assign a dual variable 7r(t;) G M (a potential) to each node v £ V, 7(a) G (a length) 
to each arc a £ E, and ^(p) G M+ to each odd fragment p £ Consider the linear program: 

minimize ^/^(vr, 7, ^) := ^ 14(0)7(0) + ^(2n(C/p) — 2)^(/)) subject to (9) 

E n 

7(a) > \/a£ E (10) 

e(/9) > ypen (11) 

7r(s) = (12) 

vr(.') = 1 (13) 

TT{v)-n{w)+j{a)+J2^ip)Xp{a)>0 ya = {v,w)eE. (14) 

n 

In fact, (|m)-()14|) is dual to linear program (^-(jlj,©. (To see this, introduce an extra arc 
(s', s), add the conservation constraints for s and s', and replace the objective (Q) by max{/(s', s)}. 
The latter generates the dual constraint 7r(s') — 7r(s) > 1. We can replace it by the equality or 
impose ^^~^^.) Therefore, 

max I/I = min '(/'(7r,7,0, (15) 

where the maximum and minimum range over the corresponding feasible solutions. 

We assert that every maximum IS-flow / achieves the maximum in ()15() . To see this, choose 
an odd barrier B = {A; Xi, . . . , X^) of minimum capacity u{B). For i = 1, . . . , A;, let C/j be the set 
of arcs from A to Xi] then pi = (Xj, Ui) is an odd fragment for G, u. Define 'k{v) to be for v £ A, 
1 for u G A', and 1/2 otherwise. Define 7(a) to be 1 for a G {A, A'), 1/2 for a G {A, M) U (M, A'), 
and otherwise, where M = - (A U ^' U Xi U . . . U Xk). Define i{pi) = 1/2 for i = 1, . . . , fe, 
and ^{p) = for the other odd fragments in (G,u). 
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One can check that ((Til) holds for all arcs a (e.g., both values ■k{w) — Tr{v) and 7(a) + 
Y^n^ip)Xpia) are equal to 1 for a = {v,w) G {A, A'), and 1/2 for a = {v,w) e iA,L) U iL,A'), 
where L ■.= V — {AU A')). Thus tt, 7, ^ are feasible. 

Using the fact that u{A, M) = u{M, A'), we observe that u • 7 = u{A, A') + u{A, M). Also 



This implies tp{Tr,^,(^) = u{B), and now the result follows from Theorem 13.51 ■ 

5 Algorithm Using a Good Pre-Solution 

Anstee |5] developed efficient methods for b-factor and b-matching problems (unweighted or 
weighted) based on the idea that a good pre-solution can easily be found by solving a corresponding 
flow problem. In this section we adapt his approach to solve the maximum IS-flow problem in a 
skew-symmetric network N = {G = (y,E),u). The algorithm that we devise is relatively simple; 
it finds a "nearly optimal" IS-flow and then makes 0(n) augmentations to obtain a maximum 
IS-flow. The algorithm consists of four stages. 

The first stage ignores the fact that N is skew-symmetric and finds an integer maximum flow 
g in N hy use of a max-flow algorithm. Then we set h{a) := {g{a) + g{a'))/2 for all arcs a £ E. 
Since div/i(s) = divg(s) /2 — divg(s') /2 = divg(s), his a maximum flow as well. Also h is symmetric 
and half-integer. Let Z be the set of arcs on which h is not integer. If Z = 0, then h is already a 
maximum IS-flow; so assume this is not the case. 

The second stage applies simple transformations of h to reduce Z. Let H = {X, Z) be the 
subgraph of G induced by Z. Obviously, for each x £ V, divh{x) is an integer, so x is incident 
to an even number of arcs in Z. Therefore, we can decompose H into simple, not necessarily 
directed, cycles Gi, . . . ,Gr which are pairwise arc-disjoint. Moreover, we can flnd, in linear time, 
a decomposition in which each cycle Gi is either self-symmetric (Gi = a{Gi)) or symmetric to 
another cycle Gj {Gi = a{Gj)). 

To do this, we start with some node vq £ X and grow in H a. simple (undirected) path 
P = {vo,ei,vi, . . . ,eg,Vq) such that the mate v'^ of each node Vi is not in P. At each step, we 
choose in H an arc e ^ eg incident to the last node Vg (e exists since H is eulerian); let x be 
the other end node of e. If none of x, x' is in P, then we add e to P. If some of x, x' is a node 
of P, Vi say, then we shorten P by removing its end part from Cj+i and delete from H the arcs 
Cj+i, . . . ,eg,e and their mates. One can see that the arcs deleted induce a self-symmetric cycle 
(when x' = Vi) or two disjoint symmetric cycles (when x = Vi). We also remove the isolated 
nodes created by the arc deletions and change the initial node vq if needed. Repeating the process 
for the new current graph H and path P, we eventually obtain the desired decomposition C, in 
0{\Z\) time. 

Next we examine the cycles in C. Each pair C, G' of symmetric cycles is canceled by sending a 
half unit of flow through G and through G', i.e., we increase (resp. decrease) h{e) by 1/2 on each 
forward (resp. backward) arc e of these cycles. The resulting function h is symmetric, and divft(x) 



Y,iMUp) - 2)e(p) = 2^Mu^) - 2) = E ^'^) - ^• 
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is preserved at each node x, whence h is again a maximum symmetric flow. Now suppose that two 
self-symmetric cycles C and D meet at a node x. Then they meet at x' as well. Concatenating 
the X to x' path in C and the x' to x path in D and concatenating the rests of C and D, we 
obtain a pair of symmetric cycles and cancel these cycles as above. These cancellations result in C 
consisting of pairwise node-disjoint self-symmetric cycles, say Ci, . . . ,Ck- The second stage takes 
0{m) time. 

The third stage transforms h into an IS-flow / whose value |/| is at most k units below \h\. 
For each i, fix a node ti in Cj and change h on Ci by sending a half unit of flow through the ij to 
t[ path in Ci and through the reverse to the t'i to ti path in it. The resulting function h is integer 
and symmetric and the divergences preserve at all nodes except for the nodes ti and t'^ where we 
have div/i(tj) = — div/i(t-) = 1 for each i (assuming, without loss of generality, that all tj's are 
different from s'). Therefore, h is, in essence, a multiterminal IS-flow with sources s,ti, . . . ,tk and 
sinks s' ,t[, . . . ,t'j^. A genuine IS-flow / from s to s' is extracted by reducing h on some /i-regular 
paths. More precisely, we add to G artificial arcs = (s, ti), i = 1, . . . , k and their mates, extend 
h by ones to these arcs and construct a symmetric decomposition P (defined in Section ^ for 
the obtained function h' in the resulting graph C (clearly h' is an IS-flow of value \h\ + k). 

Let V be the set of elementary flows in D formed by the paths or cycles which contain artificial 
arcs. Then 6 = 1 for each (P, P' , 6) £ V. Define /' := h' - Y^ix^ + X^' ■ {P, P' , 1) G ^')- Then 
/' is an IS-flow in C , and |/'| > \h'\ — 2k > \h\ — k. Moreover, since /(cj) = for i = 1, . . . ,k, 
the restriction / of /' to E is an IS-flow in C, and |/| = |/'|. Thus, |/| > \h\ — k, and now the 
facts that k < n/2 (as the nodes ti, . . . , t^, t'^, . . . , are different) and that /i is a maximum flow 
in imply that the value of / differs from the maximum IS-flow value by 0{n). The third stage 
takes 0{nm) time (the time needed to construct a symmetric decomposition of h'). 

The final, fourth, stage transforms / into a maximum IS-flow. Each iteration applies the r- 
reachability algorithm (RA) mentioned in Section [21 to the split-graph S{C~^,Uf) in order to find 
a Mj-regular s to s' path P in G"*" and then augment the current IS-flow / by the elementary flow 
(P, P' ,Suf{P)) as explained in Section |3I Thus, a maximum IS-flow in N is constructed in 0(n) 
iterations. Since the RA runs in 0{m) time (by Theorem I2.2|l . the fourth stage takes 0{nm) 
time. 

Summing up the above arguments, we conclude with the following. 

Theorem 5.1 The above algorithm finds a maximum IS-flow in N in 0{M{n,m) + nm) time, 
where M{n, m) is the running time of the max-flow procedure it applies. 

6 Shortest R- Augmenting Paths and Blocking IS-Flows 

Theorem 13 . 31 and Corollarv 13 .41 prompt an alternative method for flnding a maximum IS-flow in a 
skew-symmetric network N = {C, u), which is analogous to the method of Ford and Fulkerson for 
usual flows. It starts with the zero flow, and at each iteration, the current IS-flow / is augmented 
by an elementary flow in {C'^,Uf) (found by applying the r-reachability algorithm to S{C~^ ,Uf)). 
Since each iteration increases the value of / by at least two, a maximum IS-flow is constructed in 
pseudo-polynomial time. In general, this method is not competitive to the method of Section |SJ 
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More efficient methods involve tlie concepts of shortest r-augmenting paths and shortest 
blocking IS-flows that we now introduce. Let g be an IS-flow in a skew-symmetric network 
{H = iy,W),h). We call giW) (= '^Ze^w 9^^)) volume of g. Considering a symmetric 
decomposition D = {(Pj, P-', 5i) : i = 1, . . . ^k) oi g, we have 

9{W) = Y.{5^\P^\ +5i\Pl\:i = l,...,k)> |(7|min{|Pi| : i = 1, . . . , fc}. 

This implies 

g{W)>\g\T-d\sis{H,h){s.s'), (16) 

where r-distj7'(x, y) denotes the minimum length of a regular x to y path in a skew-symmetric 
graph H' (the regular distance from x to y). We say that an IS-flow g is 

(i) shortest if H16() holds with equality, i.e., some (equivalently, any) symmetric decomposition 
of g consists of shortest /i-regular paths from s to s'; 

(ii) totally blocking if there is no {h — (7)-regular path from s to s' in H, i.e., we cannot augment 
g using only residual capacities in H itself; 

(iii) shortest blocking if g is shortest (as in (i)) and 

T^-dists^H,h-g){s, s') > r-dist5(^^_;,)(s, s'). (17) 

Note that a shortest blocking IS-flow is not necessarily totally blocking, and vice versa. 

Given a skew-symmetric network N = (G, u), the shortest r-augmenting path method (SAPM), 
analogous to the method of Edmonds and Karp [7] for usual flows, starts with the zero flow, and 
each iteration augments the current IS-flow / by a shortest elementary flow g = {P, P' ,Su^.(P)). 

The shortest blocking IS-flow method (SBFM), analogous to Dinits' method starts with 
the zero flow, and each phase (big iteration) augments the current IS-flow / by performing the 
following two steps: 

(PI) Find a shortest blocking IS-flow g in {G^,Uf). 

(P2) Update f ■= f ® g. 

Both methods terminate when / no longer admits r-augmenting paths (i.e., g becomes the 
zero flow). The following observation is crucial for our methods. 

Lemma 6.1 Let g be a shortest IS-flow in {G~^,Uf), and let f ■= f (B g- Let k and k' be the 
minimum lengths of r-augmenting paths for f and f , respectively. Then k' > k. Moreover, if g 
is a shortest blocking IS-flow, then k' > k. 

Proof. Take a shortest u//-regular path P from s to s' in G"*". Then |P| = k' and g' = {P, P' , 1) 
is an elementary flow in {G'^^Uf). 

Note that supp(5f) does not contain opposed arcs a = {x,y) and b = {y,x). Otherwise 
decreasing g by one on each of a,b,a',b' (which are, obviously, distinct), we would obtain the 
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IS-flow g in {G~^,Uf) such that \g\ = \g\ and g{E~^) < g{E'^), which is impossible because 
9{E~^) ^ ^1^1 and g{E^) = k\g\. This imphes that each arc a in the set Z := {a E : g{a^) = 0} 
satisfies 

Uf'{a) = Uf{a)-g{a). (18) 

If supp(g') C Z, then g' is a feasible IS-flow in (G+,m/) (by whence k' = g'{E+)/\g'\ > k. 

Moreover, if, in addition, g is a shortest blocking IS-flow, then ((T7|l and the fact that g' <Uf — g 
(by (dHI)) imply k' > k. 

Now suppose there is an arc e G E~^ such that g'{e) > and g{e^) > 0. For each a £ i?"*", put 
A(a) := max{0,5(a) + g''(a) — g{a^) — g'{a^)}. One can check that A(a) < ?ij(a) for all arcs a and 
that div;^(u) = for all nodes v ^ s, s'. Therefore, A is an IS-flow in {G~^ ,Uf) with |A| = \g\ + \g'\ = 
\g\ + 2. Also X{E+) < g{E+) + g'{E+) since for the e above, A(e) + A(e^) < g'{e) + g{e^). We 
have 

2k' = g'{E+) > X{E+) - g{E+) > k{\g\ + 2) - k\g\ = 2k, 

yielding k' > k. ■ 

Thus, each iteration of SAPM does not decrease the minimum length of an r-augmenting 
path, and each phase of SBFM increases this length. This gives upper bounds on the numbers of 
iterations. 

Corollary 6.2 SAPM terminates in at most (n — l)m iterations. 

(This follows by observing, in the proof of Lemma l6.ll that on the iterations with the same 
length of shortest r-augmenting paths, the subgraph of induced by the arcs contained in such 
paths is monotone nonincreasing, and each iteration reduces the capacity of some arc of this 
subgraph, as well as the capacity of its mate, to zero or one.) 

Corollary 6.3 SBFM terminates in at most n — 1 phases. 

As mentioned above, SBFM can be considered as a skew-symmetric analog of Dinits' blocking 
flow algorithm. Recall that each phase of that algorithm constructs a blocking flow in the sub- 
network H formed by the nodes and arcs of shortest augmenting paths. Such a network is acyclic 
(moreover, layered), and a blocking flow in H is easily constructed in 0{nm) time. 

The problem of finding a shortest blocking IS-flow ((PI) above) is more complicated. Let H 
be the subgraph of G^ formed by the nodes and arcs contained in shortest uj-regular s to s' 
paths. Such an H need not be acyclic (a counterexample is not difficult). In Section El we will 
show that problem (PI) can be reduced to a seemingly easier task, namely, to finding a totally 
blocking IS-flow in a certain acyclic network (H,h). Such a network arises when the shortest 
r-path algorithm from 16 is applied to the split-graph S{G~^,Uf) with unit arc lengths. First, 
however, we need to tell more about the r-reachability and shortest r-path algorithms from |16j . 

7 Properties of Regular and Shortest Regular Path Algorithms 

In this section we exhibit certain properties of the algorithms in (21) referring the reader to that 
paper for details. We also establish an additional fact fLemma l7.4() . which will be used later. 
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Figure 3: Fragment trimming example 



7.1 The Regular Reachability Algorithm (RA) 

Let F = {V, E) be a skew-symmetric graph with source s and sink s' = (j{s) (as before, a is the 
symmetry map). A fragment (or an s-fragment) in F is a pair (j) = (V^,e0 = {v,w)), where is 
a self-symmetric set of nodes of F with s and is an arc entering V^, i.e., v ^ B w (cf. 
the definition of odd fragments in Section llj. We refer to e,^ and as the base and anti-base arcs 
of (f), respectively. Let us say that the fragment is well-reachable if 

(i) for each node x S V^, there is an r-path from u) to x in the subgraph induced by (and 
therefore, an r-path from x to w' = a{w)), and 

(ii) there is an r-path from s to w disjoint from Vr- 

The trimming operation applied to (p (which is analogous to shrinking a blossom in matching 
algorithms) transforms F by removing the nodes of — {w, w'} and modifying the arcs as follows. 

(Tl) Each arc a = {x,y) £ E such that either x,y €z V — or a = e^j) or a = e'^ remains an arc 
from X to y. 

(T2) Each arc {x, y) £ E — {e^} that leaves is replaced by an arc from w to y, and each arc 
(x, y) G E — {e^} that enters is replaced by an arc from x to w'. 

(T3) Each arc with both ends in is replaced by an arc from w to w' . 

(A variant of trimming deletes all arcs in (T3).) The image of an arc a in the new graph is denoted 
again by a (so its end nodes can be changed, but not its name). Figure 01 illustrates fragment 
trimming. The new F is again skew-symmetric. 

The algorithm RA relies on the following property. 

Statement 7.1 |16j If (j) is a well-reachable fragment, then trimming (j) preserves the existence 
(or non-existence) of a regular path from s to s' . 
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RA searches for a regular s to s' path in T, starting with the trivial path P = s. Each iteration 
either increases the current r-path P, or reveals a well-reachable fragment and trims it, producing 
the new current graph F and accordingly updating P. It terminates when either an s to s' r-path 
P or an s-barrier B in the final graph T is found (cf. Theorem I2.1j) . The postprocessing stage 
extends P into a regular s to s' path P of the initial T (cf. Statement 17. 1|) in the former case (the 
path restoration procedure) and extends B into a barrier B of the initial T in the latter case (the 
barrier restoration procedure). 

The fragments of current graphs revealed by RA determine fragments of the initial F in a 
natural way; all fragments are well-reachable. Moreover, the set <I> of these fragments of the 
initial F is well-nested. This means that 

(Fl) for distinct (t>,if; either V^/, C or C or V^i, n = $ , and 

(F2) for G ^, if K/) C and G S{V(j)) then = e^, and if n = and G '5(K/>) 
then 60 6{V^). 

(Recall that for X <^V, S{X) is the set of arcs with one end in X and the other in y — X.) Let 
us say that a path i? in F is compatible with $ if for each (j) £ ^, p := |i?n (5(V^)| < 2, and if p = 2 
then R contains exactly one of e^f,, e'^. The following additional properties (relying on (T1)-(T3)) 
are important: 

(19) any regular s to s' path P in the final graph F is extendable to a regular s to s' path 
P compatible with $ in the initial F, and the path restoration procedure applied to P 
constructs such a P in 0(|P| -|- d) time, where d is the total size of maximal fragments in <I> 
traversed by P\ 

(20) for each G <I> and each arc a ^ e'^ leaving l^, there exists a compatible with $ r-path 
Q(f){a) with the first arc e<^, the last arc a and all intermediate nodes in V^; such a path can 
be constructed by (a phase of) the path restoration procedure in 0(1^^1) time. 

A fast implementation of RA (supported by the disjoint set union data structure of JH]) runs 
in linear time, as indicated in Theorem 12.21 

7.2 The Shortest Regular Path Algorithm (SRA) 

We now consider the shortest regular path problem (SRP) in a skew-symmetric graph F = (y, E) 
with nonnegative symmetric lengths i{e) of the arcs e £ E: find a minimum length regular path 
from s to s' . One may assume that s' is r-reachable from s. The dual problem involves above- 
mentioned fragments. Define the characteristic function X(t> of ^ fragment (j) = (V^, e,^) by 




(21) 
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(Compare with For a function ir : V ^M. (of node potentials) and a nonnegative function ^ 
on a set $ of fragments, define the reduced length of an arc e = (x, y) to be 



iiie) := i{e) + 7r(x) - ^(y) + J] e(</')x^(e). 



An optimality criterion for SRP can be formulated as follows. 

Theorem 7.2 jJB] yl regular path P from s to s' is a shortest r-path if and only if there exist a 
potential tt : 1/ — > M, a set ^ of fragments, and a positive function £^ on ^ such that 



The shortest r-path algorithm (SRA) from JB] implicitly maintains tt, in the input graph 
r and iteratively modifies the graph by trimming certain fragments. Let be the subgraph of 
the current F with the same set of nodes and with the arcs having zero reduced length, called the 
current 0-subgraph (recall that the arcs of the current graph are identified with the corresponding 
arcs of the initial one). Each iteration applies the above r-reachability algorithm RA to search for 
a regular s to s' path in F*^. If such a path is found, the algorithm terminates and outputs this path 
to a postprocessing stage. If such a path does not exist, then, using the barrier B = {A; Xi, . . . , Xk) 
in F'' constructed by RA, the iteration trims the fragments determined by the sets Xi and updates 
IT, <I?,^, modifying F*^. The reduced lengths of the arcs within the newly and previously extracted 
fragments, as well as of their base and anti-base arcs, are not changed. 

Let F and F denote the final graph F and the 0-subgraph in it, respectively, and P the regular 
s to s' path in F found by the algorithm. Let F^ stand for the 0-subgraph of the initial graph F 
(concerning the reduced arc lengths determined by the resulting vr,^,^). We call F'^ and F the 
full and trimmed 0-graphs, respectively. The postprocessing stage applies the path restoration 
procedure of RA to extend P into a regular s to s' path P in F*^, in time indicated in (fT9|) . It also 
explicitly constructs F*^ (in linear time). 

Note that any s to s' r-path or r-cycle Q in F compatible with $ satisfies x'^ ' = for each 
G <I>. By ()19() . any s to s' r-path P in f'' is extendable to an s to s' r-path P in F*^ compatible 
with <I>. Therefore, P is shortest, by Theorem 17. 21 

Theorem 7.3 ^] For nonnegative symmetric arc lengths i, SRA runs in 0(m- log n) time, and 
in 0{m^/\og L) time when i is integer-valued and L is the maximum arc length. Furthermore, 
the algorithm constructs (implicitly) 7r,<I>,^ as in Theorem \ 7. S[ where tt is anti-symmetric (i.e., 
-7r(x) = —Tr{x') for all x ^ V), and constructs (explicitly) the trimmed 0-graph F and the full 
0-graph such that: 

(Al) <I> is well-nested (obeys (F1)-(F2)J and consists of well-reachable fragments in T^; in par- 
ticular, ^1(6^) = for each (/> G <I>; 



iiie) > for each e e E; 
£|(e) = for each e e P; 



(22) 
(23) 
(24) 



• X0 = foi^ each (/) G <I>. 
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(A2) $ satisfies liiy\) and with r'^,T instead ofT,T; in particular, any regular s to s' path 
ofT^ is (efficiently) extendable to a shortest regular s to s' path in {T,£). 

(Note that the anti-symmetry of vr and the symmetry of i and x<j) all G <^ imply that 
the reduced length function £| is symmetric. Therefore, the graphs and are indeed skew- 
symmetric.) Let <I)™^^ denote the set of maximal fragments in <I>. The sets for (p £ (j)max g^^.^ 
pairwise disjoint (by (Fl)), and the graph can be directly obtained from by simultaneously 
trimming the fragments in (J)™^^. 

In the next section we will take advantage of the relationship between r-paths in T and 

shortest r-paths in {T,£) indicated in (A2). Another important property of T is as follows. 

Lemma 7.4 // the length £{C) of every cycle C inV is positive, then is acyclic. In particular, 
r is acyclic if all arc lengths are positive. 

Proof. Suppose contains a (not necessarily regular) simple cycle C. In view of (|2Uj) . C is 
extendable to a cycle C of compatible with $. Then x*^ ' Xci> = ^ for all </) € <I*. This implies 
that the original length £{C) and the reduced length ^I(C) are the same (since the changes in 
£| due to vr cancel out as we go around the cycle). Since all arcs of C have zero reduced length, 
£(C) = ^I(C) = 0. This contradicts the hypotheses of the lemma. ■ 



8 Reduction to an Acyclic Network and Special Cases 

We continue the description of the shortest blocking IS-flow method (SBFM) for solving the 
maximum IS-flow problem in a network N = {G = {V,E),u) begun in Section |HJ Let / be a 
current IS-flow in A^. We show that the task of finding a shortest blocking IS-flow g in {G^ ,Uf) 
(step (PI) of a phase of SBFM) reduces to finding a totally blocking IS-flow in an acyclic network. 

Build the split-graph F = S{G'^,Uf) and apply the above shortest regular path algorithm to 
F with the all-unit length function £ on the arcs. It constructs </), as in Theorems 17.21 and 17.31 
taking 0(m) time (since L = 1). SRA also constructs the trimmed 0-graph F , the main object 
we will deal with. By Lemma 177^ f'^ is acyclic. Also the following property takes place. 

Lemma 8.1 Let a E be an arc with Uf{a) > 1, and let ai,a2 be the corresponding split-arcs 
in F. Then £t{ai) = £^{02). Moreover, none of 01,02 can be the base or anti-base arc of any 
fragment in <I>. 

Proof. Since 01,02 are parallel arcs, for each G oi enters (resp. leaves) if and only if 
02 enters (resp. leaves) V^. This implies that ^l(ai) / ^I(a2) can happen only if one of 01,02 is 
the base or anti-base arc of some fragment in <I>. Suppose oi £ {e^, e'^} for some (j) £ ^ (the case 
^2 G is similar). Then ^l(ai) = (by (Al) in Theorem 17. 3|) . Using property (F2) from 

Section [7| (valid as <I> is well- nested), one can see that 02 is not the base or anti-base arc of any 
fragment in $. Therefore, Xi>{(^2) < Xip{(^i) ^^v all -0 G yielding it (02) < £|(ai). Moreover, 
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the latter inequality is strict because x</>(fl2) = — 1 < 1 = Xtpio-i) and ^((/)) > 0. Now zi{ai) = 
implies £|(a2) < 0, contradicting (221) • ■ 

Let C be the set of (images of) zero reduced length arcs of T. Lemma 18.11 implies 
that the base arc of each fragment (/> G <I> in F is generated by an arc e G E^ with Uf{e) = 1. 
We can identify these e and and consider (p as a fragment of as well. One can see that 
r° is precisely the split-graph for {H,h), where H = {V,E^) is obtained from H = {V,E^) by 
trimming the maximal fragments in <I>, and h is the restriction of Uf to E . 

Based on the property of each fragment to have unit capacity of the base arc, we reduce step 
(PI) to the desired problem, namely: 

(B) Find a totally blocking IS- flow in {H, h). 

To explain the reduction, suppose we have found a solution 'g to (B). For each maximal 
fragment (/> in $ with € supp(^), we have ^(e,^) = 1; therefore, exactly one unit of flow goes 
out of the head of e^, through an arc a G ^ say. We choose the path Q = Q^{a) as in (|20j) to 
connect and a in (the subgraph on V(f) of) H and then push a unit of flow through Q and a 
unit of flow through the symmetric path Q'. Doing so for all maximal fragments (j), we extend g 
to an IS-flow g in {H, h), where h is the restriction of to E^. Moreover, g is a shortest blocking 
IS-flow in (G+,n/). 

Indeed, the fact that the chosen paths Q have zero reduced length and are compatible with 
$ implies that a symmetric decomposition of g consists of shortest ttj-regular paths (cf. (A2) in 
Theorem 17. so g is shortest. Also cannot contain a (n/ — g)-iegulaT s to s' path R of length 
g{E^)/\g\. For such an R would be a path in H compatible with <I> (in view of Theorem 17. 2 j) : 
then the arcs of R occurring in H should form an (h — ^)-regular s to s' path in it, contrary to 
the fact that g is totally blocking. 

Since each path Q^{a) is constructed in OdV^I) time (by (jSHl)), and the sets of maximal 
fragments (j) are pairwise disjoint, the reduction to (B) takes linear time. 

Lemma 8.2 A totally blocking IS-flow in (H, h) can be extended to a shortest blocking IS-flow 
in {G~^,Uf), in 0[m) time. ■ 

Corollary 8.3 SBFM solves the maximum IS-flow problem in 0{qT{n,m) + qm) time, where q 
is the number of phases (q < n) and T{n,m) is the time needed to find a totally blocking IS-flow 
in an acyclic network with at most n nodes and m arcs. 

Clearly T{n,m) is O(m^), as a totally blocking flow can be constructed by 0{m) applications 
of the regular reachability algorithm; this is slower compared with the phase time 0{nm) in 
Dinits' algorithm. However, we shall show in the next section that problem (B) can be solved in 
0{nm) time as well. Moreover, the bound will be better for important special cases. 

Next we estimate the number of phases. For the standard max-flow problem, the number of 
phases of Dinits' algorithm becomes significantly less than n in the cases of unit arc capacities 
and unit "node capacities" . To combine these into one case, given a network N = [G = {V, E),u) 
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with integer capacities u, for a node x € V, define the transit capacity u(x) to be the minimum 
of values Ey:ix,y)&Eu{x,y) and J2y:{y,x)eEu{y,x). Define 

A := A{N) := ^{u{x) :xeV-{s, s'}). 

As shown in [Hll^), the number q of phases of the blocking flow method does not exceed 2\/A. In 
particular, if m = 1 then q = 0{y/m), and if the transit capacities u{x) of all nodes x ^ s,s' (inner 
nodes) are ones, e.g., in the case arising from the bipartite matching problem, then q = 0{^/n). 
A similar argument works for skew-symmetric networks (see also for a special case). 

Lemma 8.4 The number of phases of SBFM is at most min{n, 2\/A}. 

Proof. After performing d := \/A phases, the r-distance from s to s' in the network A^' = (G"*", Uf) 
for the current IS-flow / becomes greater than d, by Lemma Id. II Let /* be a maximum IS-flow 
in A^, and let g be defined as in the proof of Theorem 13.31 Then g is a feasible IS-flow in N' and 
Iffl — \ f*\ ~ I/I- We assert that \g\ < d, which immediately implies that the number of remaining 
phases is at most d/2, thus proving the lemma. To see this, take a symmetric decomposition T> of 
g consisting of elementary flows (P, P' , 6) with 6 = 1. Let T>' be the family of s to s' paths P, P' 
in P; then > \g\. It is easy to see that Uf{x) = u{x) for each inner node x. Each path in T>' 
contains at least d inner nodes, and therefore, it uses at least d units of the total transit capacity 
of inner nodes of N'. So we have dl^^'l < A(A). This implies 1^1 < d. ■ 

9 Finding a Totally Blocking IS-Flow in an Acyclic Network 

Our aim is to show the following. 

Theorem 9.1 A totally blocking IS-flow in an acyclic skew- symmetric graph with 0{n) nodes, 
0{m) arcs, and unit arc capacities can be found in 0{n-\- m) time. 

This together with Corollary 18.31 and Lemma 18.41 yields the following result for the shortest 
blocking IS-flow method. 

Corollary 9.2 In case of a network N with unit arc capacities, SBFM can be implemented so 
that it finds a maximum IS-flow in 0{m^ A(A)) time (assuming n = 0{m)). In particular, if 
the indegree or outdegree of each node is at most one, then the running time becomes 0{^/nm). 

In the second half of this section we will extend Theorem 19. II to general capacities, in which 
case the phase time will turn into 0{nm), similarly to Dinits' algorithm. 

For convenience we keep the original notation for the network in question. Let G = (V, E) be 
a skew-symmetric acyclic graph with source s and the capacity u{e) = 1 of each arc e £ E. One 
may assume that each node belongs to a path from s to cr(s). 

First of all we make a reduction to the maximal balanced path-set problem (MBP) stated in 
the Introduction. Since G is acyclic, one can assign, in linear time, a potential function vr : y — > Z 
which is antisymmetric (vr(x) = — -7r(f7(x)) for each x ^V) and increasing on the arcs (7r(y) > 7r(x) 
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for each (x, y) £ E). (Indeed, a function g : F — > Z increasing on the arcs is constructed, in hnear 
time, by use of the standard topological sorting. Now set Tr{v) := qiv)—q{a{v)), v G V.) Subdivide 
each arc (x, y) with 7r(x) < and 7r(y) > into two arcs (x, z) and (z, y) and assign zero potential 
to z. The new graph G, with 0{m) nodes and 0{m) arcs, is again skew-symmetric, and the 
problem remains essentially the same. 

Let r be the subgraph of the new G induced by the nodes with nonnegative potentials. Then 
r U o"(r) = G and T n o"(r) = (Z, 0), where Z is the self-symmetric set of zero potential nodes of 
G. Also r contains (t{s). 

Clearly every s to (t{s) path P of G meets Z at exactly one node z, which subdivides P into 
an s to z path R' in o"(r) and a z to a{s) path Q in T. Then P is regular if and only if <j{R') 
and Q are arc-disjoint. Conversely, let Q,R be two arc-disjoint Z to (j[s) paths in F beginning 
at symmetric nodes in Z. Then the concatenation (t{Q) ■ R (as well as (j{R) • Q) is a regular s to 
a{s) path of G if and only if Q and R are arc-disjoint. 

This shows that our particular totally blocking IS-flow problem is reduced, in linear time, 
to MBP with T,a(s),Z (in fact, the problems are equivalent). Theorem 19.11 is implied by the 
following. 

Theorem 9.3 MBP is solvable in linear time. 

We devise an algorithm for MBP and prove Theorem 19.31 Let the input of MBP consist of 
an acyclic graph T = {X, U), a sink t and a source set Z with a map (involution) a : Z ^ Z 
giving a partition of Z into pairs. We say that two arc-disjoint Z to t paths Q,R beginning at 
"symmetric" sources z,a{z) form a good pair, and say that a collection of pairwise arc-disjoint Z 
to s' paths in r is a balanced path-set if its members can be partitioned into good pairs. So the 
task is to find a maximal (or "blocking") balanced path-set. 

Each iteration of the algorithm will reduce the arc set of F and, possibly, the set Z, and we 
sometimes will use index i for objects in the input of i-th iteration. So Fi = (Xi, Ui) is the initial 
graph. Without loss of generality, one may assume that initially each source has zero indegree 
and 

(CI) each node of F lies on a path from Z to t, 

and will maintain these properties during the algorithm. 

The iteration input will include a path D from a certain node of F to t, called the pre-path. 
Initially, D is trivial: D = t. The nodes of F not in Z U {t} are called inner. The current F may 
contain special inner nodes, called complex ones. They arise when the algorithm shrinks certain 
subgraphs of F; the initial graph has no complex nodes. The adjacency structure of F is given 
by double-linked lists Ix and Ox of the incoming and outgoing arcs, respectively, for each node x. 
The arc set of a path P is denoted by E{P). 

An i-th iteration begins with extending D to a Z to t path P in a natural way; this takes 
0(|P| — \D\) time. Let z be the first node of P. Then we try to obtain a good pair by constructing 
a path from z' = cr(z) to t, possibly rearranging P. By standard arguments, a good pair for z, z' 
exists if and only if there exists a path A from z' to t, with possible backward arcs, in which the 
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forward arcs belong to U —E{P) and the backward arcs belong to E{P), called an augmenting path 
w.r.t. P. For certain reasons, we admit A to be self-intersecting in nodes (but not in arcs). Once 
A is found, the symmetric difference E{P)AE{A) gives a good pair Q,R (taking into account 
that r is acyclic). 

To search for an augmenting path, we replace each arc e = (x, y) £ E{P) by the reverse arc 
e = {y,x); let T = {X,U) be the resulting graph, and P the t to Z path reverse to P. Thus, we 
have to construct a (directed) path from z' to t in F or establish that it does not exist. 

To achieve the desired time bound, we apply a variant of depth first search which we call here 
transit depth first search (TDFS) (such a search procedure was applied in JH])- The difference 
from the standard depth first search (DFS) is as follows. When scanning a new outgoing arc {x, y) 
in the list Ox of a current node x, if y has already been visited, then DFS stays at x. In contrast, 
TDFS moves from x to y, making y the new current node. Both procedures maintain the stack of 
arcs traversed only in forward direction and ordered by the time of their traversal. If all outgoing 
arcs of the current node x are already scanned, then the last arc {w, x) of the stack is traversed 
in backward direction and w becomes the new current node. We refer to the path determined by 
the stack, from the initial node to the current one, as the active path. Note that in case of TDFS 
the active path may be self- intersecting (while it is simple in DFS). 

We impose the condition that the outgoing arc lists of F are arranged so that 

(C2) for each node x ^ z oi P, the arc e of P leaving x is the last element of Ox- 

This guarantees that TDFS would scan e after all other outgoing arcs of x (i.e., the arcs of P are 
ignored as long as possible). 

At an iteration, we apply TDFS to F starting from z' as above. The search terminates when 
either it reaches t or it returns to z' having all arcs of Oz' traversed. In the first case (breakthrough) 
the final active path ^ in F determines the desired augmenting path ^ in F, and we create a good 
pair Q, R as described above. In the second case, the non-existence of a good pair for the given 
z, z' is declared. Consider both cases. 

Breakthrough case. Delete from F the arcs of Q,R and then delete all the nodes and arcs 
that are no longer contained in Z to t paths (thus maintaining (CI)). This is carried out by 
an obvious cleaning procedure in 0{q) time, where q is the number of arcs deleted. If Q or i? 
contains a complex node, the iteration finishes by transforming Q, R into a good pair of paths of 
the initial graph; this is carried out by a path expansion procedure which will be described later. 
The obtained F, Z form the input of the next iteration, and the new pre-path D is assigned to be 
the trivial path t. If F vanishes, the algorithm terminates. The following observation is crucial 
for estimating the time bound. 

Lemma 9.4 Let q be the number of arcs deleted at an iteration with a breakthrough. Then, 
excluding the path expansion procedure if applied, the iteration runs in 0{q) time. 

Proof. Let W be the set of arcs of F traversed by TDFS on the iteration, and W the corresponding 
set in F, i.e., W = {e € U : e G W or e € W}. The iteration runs in 0{q + \W\) time, taking into 
account that each arc of P not in Q U i? is contained in W. Therefore, it suffices to show that no 
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arc from W remains in the new graph F. Suppose this is not so. Then there is a Z to t path L of 
the old r that meets W but not E{Q) U E{R) (as the arcs of Q U R are deleted). Let e = {x, y) 
be the last arc of L in W. Let b = {y, w) be the next arc of L (it exists since y = t would imply 
that e is in A but not in P, whence e belongs to QU R). Then b ^WU E{Q) U E{R), by the 
choice of e. Two cases are possible. 

(i) e is in P. Then e = {y,x) £ W. According to condition (C2), at the time TDFS traversed 
e from y to x all arcs of T leaving y had already been traversed. So b is not in F, implying 
b £ E{P) — W. Then 6 is in Q U -R; a contradiction. 

(ii) e is not in P. Then e £ W and e does not lie on the final active path A (otherwise e is in 
QUR). Therefore, TDFS traversed e in both directions. To the time of traversal of e in backward 
direction, from y to x, all arcs of F leaving y have been traversed (at this point the difference 
between TDFS and DFS is important). So b is not in F, whence b £ E{P). Now b ^ W implies 
that 6 is in Q U i?; a contradiction. ■ 

Non-breakthrough case. Let Y be the set of nodes visited by TDFS. Then no arc of F leaves 
Y. Therefore, in view of (CI), 

(25) the set of arcs of F leaving Y consists of a unique arc a = {v, w), this arc lies on P, and 
the nodes of the part of P from z to v are contained in Y. 

Since no arc of F enters Z, we also have 

YnZ = {z,z'}. (26) 

We reduce F by shrinking its subgraph Fy = (Y, Uy ) induced by Y into one node; the formed 
complex node vy is identified with v. We call v the root of Fy and store Fy. The list of arcs 
entering vy in the new graph is produced by simply merging the lists Ix for x £ Y from which 
the arcs occurring in Fy are explicitly removed, using the lists Oy for y £Y . (We do not need 
to correct the outgoing arc lists Ox for x y explicitly, as we explain later.) Thus, to update F 
takes time linear in \Uy\. By H25() . a is the only arc leaving vy in F. 

From (CI) and ^ it follows that 

(27) the new graph F is again acyclic, and for each x £Y , there is a path Py{x) from x to 
the root v in Fy. 

In view of ()26|) . the set Z is updated as Z := Z — {z, z'}. If there is at least one arc entering 
Vy, then the new graph F and set Z satisfy (CI) and form the input of the next iteration. The 
new pre-path D is assigned to be the part of P from the formed complex node vy to t. If no arc 
enters vy, we finish the current iteration by removing the nodes and the arcs not contained in Z 
to t paths. This further reduce F and may reduce Z and shorten D. 

One can see that the set Uy is exactly W. This and the construction of pre-paths imply the 
following. 

Lemma 9.5 Let q be the number of arcs deleted by an i-th iteration without a breakthrough. Then 
the iteration runs in 0{q + max{0, — |-Dj|}) time. ■ 
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As mentioned above, we do not need to explicitly correct the outgoing arc lists Ox for x 
(this would be expensive). Let V be the current set of all complex nodes created from the beginning 
of the algorithm. We take advantage of the following facts. First, the elements of V that are nodes 
of the current graph (the maximal complex nodes) lie on the current pre-path D. Second, at an 
iteration with a breakthrough, all complex nodes are removed. Third, at an iteration without a 
breakthrough, the subgraph Ty forming the new complex node vy contains a subpath of P from 
its beginning node (by ()25() '). and the cleaning procedure (if applied at the iteration) deletes a 
part of the updated P from its beginning node as well. Therefore, one can store V as a tree in a 
natural way and use the disjoint set union data structure from to maintain V. This enables 
us to efficiently access the head vy of any arc e = {x,vy) when e is traversed by TDFS (with 
0(1) amortized time per one arc). 

To complete the algorithm description, it remains to explain the path expansion procedure to 
be applied when an iteration with a breakthrough finds paths Q, R contained complex nodes. 
It proceeds in a natural way by recursively expanding complex nodes occurring in the current 
Q,i? into the corresponding paths Py{x) as in H27() and building Py{x) into Q or R (this takes 
0{\Py{x)\) time). The arc sets of subgraphs Ty extracted during the algorithm are pairwise 
disjoint, so the total time for all applications of the procedure is 0{m). 

Thus, we can conclude from Lemmas 19.41 and 19.51 that the algorithm runs in 0{m) time, 
yielding Theorem 19. HI 

In the rest of this section we extend the above approach and algorithm (Algorithm 1) to a 
general case of acyclic {G,u). The auxiliary graph F = (X, U) and the set Z are constructed as 
above, and the capacity u{e) of each arc e G U is defined in a natural way. We call an integer Z 
to t flow g in (F,u) balanced if the flow values out of "symmetric" sources are equal, i.e., 

divg(z) = divg{a{z)) for each z G Z, 

and blocking balanced if there exists no balanced flow g' satisfying g ^ g' > g (taking into account 
that F is acyclic). Then the problem of finding a totally blocking IS-flow in (G, n) is reduced to 
problem BBF: find a balanced blocking flow for F, u, Z, t. 

Algorithm 2 will find a balanced blocking flow g in the form g = aix^^ +aix^^ + • • ■ + arX^'^ + 
OrX^^'j where each ai is a positive integer, Qi is a path from some z G Z to t, and Ri is a path from 
cr{z) to t. It iteratively constructs pairs Qi,Ri for current T,u,Z, assigns the weight to them 
as large as possible, and accordingly reduces the current capacities as n := n — Oix'^'' — OiiX^^- All 
arc capacities in F are positive: once the capacity of an arc becomes zero, this arc is immediately 
deleted from F. 

Each pair Qi, Ri is constructed as in Algorithm 1 when it is applied to the corresponding split- 
graph S = S{T,u). More precisely (cf. Section S is formed by replacing each arc e = {x,y) of 
F by two parallel arcs (split-mates) ei,e2 from x to y with the capacities \u{e)/2] and [u{e)/2\, 
respectively. When u{e) = 1, 62 vanishes in S, and ei is called critical. The algorithm maintains 
S explicitly. The desired pair Qi, Ri in (F, u) is determined by a good pair in in a natural way. 

The main part of an iteration of Algorithm 2 is a slight modification of an iteration of Algo- 
rithm 1. The difference is the following. While Algorithm 1 deletes all arcs of the paths Q,R 
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found at an iteration, Algorithm 2 deletes only a nonempty subset B of arcs \n Q\J R (concerning 
the graph S) including all critical arcs in these paths. One may think that Algorithm 2 essentially 
treats with a graph S (ignoring (r,u)) in which some disjoint pairs of parallel arcs (analogs of 
split-mates) are distinguished and the other arcs are regarded as critical, and at each iteration, 
the corresponding subset B C E{Q) U E{R) to be deleted is given by an oracle. Emphasize that 
the unique arc leaving a complex node is always critical. Therefore, each complex node m QU R 
will be automatically removed. Computing Oj's and other operations of the algorithm beyond the 
work with the graph S do not affect the asymptotic time bound. 

We now estimate the complexity of an iteration of Algorithm 2. In case without a break- 
through, properties H25() .(|26 p .()27 |) and Lemma 19.51 fwith S instead of F) remain valid. Note that 
the arc a in H25() is critical (since it is a unique arc leaving Y); therefore, the arc leaving the cre- 
ated complex node is critical. Our analysis of the breakthrough case involves the subset W2 C W 
of arcs traversed by TDFS in both directions (where W is the set of all traversed arcs in the 
corresponding auxiliary graph S). Let W2 be the corresponding set in S. 

Lemma 9.6 Suppose an iteration of Algorithm 2 results in a breakthrough. Let e = {x, y) be an 
arc of S such that e E W2 or e £ E{P) n W . Then any x to t path L in S starting with the arc e 
contains a critical arc in Q U R (and therefore, e vanishes in the new graph S). 

Proof. Suppose this is not so and consider a counterexample (e,L) with \L\ minimum. Let Aq 
be the active path in S just before the traversal of e or e from y to x, and ^0 the corresponding 
(undirected) path in S. At that time, for the set Oy of arcs of S leaving y, 

(28) all arcs in Oy except e (in case e E E(P)) are already traversed 

(in view of condition (C2)). Let b = (y, w) be the second arc of L (existing as y = t is impossible). 
By ((2HI), if b is not in P, then b £ W. Also b ^ W2 and b E{P) n W (otherwise the part of 
L from y to t would give a smaller counterexample). This implies that b belongs to Q U R, and 
therefore, b is not critical. Let b' be the split-mate of b. Considering the path starting with b' and 
then following L from w to t (which is smaller than L), we similarly conclude that b' W2 and 
b' ^ E{P) n W . To come to a contradiction, we proceed as follows. 

The fact that S is acyclic implies that the symmetric difference (on the arcs) of P and Aq is 
decomposed into a path from Z to t and a path from Z to y; therefore, E{P)AE{Aq) contains 
at most one arc a leaving y. This and imply that all arcs in Oy Pi Oy except, possibly, a 
have been traversed twice; so they are in W2- Hence, one of 6, b' must be in P; let for definiteness 
b E E{P) (then b' E{P)). 

Now b^W implies b E E{P) - E{Ao), and b' e W - W2 implies b' E £^(^0) - E{P). Thus, 
both arcs 6,6' leaving y are in E{P)AE{Aq); a contradiction. ■ 

The running time of an iteration with a breakthrough is 0(1^*1 + \W\ + q), where q is the 
number of arcs deleted from S. Lemma [9.61 allows us to refine this bound as 0(|(5| -|- -|- q)- 
Combining this with Lemma 19.51 we can conclude that, up to a constant factor, the total time of 
Algorithm 2 is bounded from above by m plus the sum E of lengths of paths Qi, Ri, . . . ,Qr, Rr 
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in the representation of the flow g constructed by the algorithm. Since |Qj|, < ?t. and r < 2m 
(as each iteration decreases the arc set of 5*), S is 0{nm). Also S does not exceed the sum of 
the transit capacities u{x) of inner nodes x of F (assuming, without loss of generality, that no arc 
goes from s to s'). Thus, Theorem I9.1l is generalized as follows. 

Theorem 9.7 For an acyclic capacitated skew- symmetric network N with 0{n) nodes and 0{m) 
arcs, a totally blocking IS-flow can be found in 0(min{m + A{N),nm}) time. 

Together with Corollary 18.31 and Lemma 18.41 this yields the desired generalization. 

Corollary 9.8 SBFM can be implemented so that it finds a maximum IS-flow in an arbitrary 
skew- symmetric network N in 0(min{n^m, yK{N){m + A(A^))}) time. 

10 Applications to Matchings 

Apply the reduction of the maximum u-capacitated b-matching problem (CBMP) in a graph 
G' = {V',E') to the maximum IS-flow problem in a network N = (G = {V, E), a,u, s); see 
Section 121 The best time bound for a general case of CBMP is attained by applying the algorithm 
of Section |SJ Theorem 15.11 implies the following. 

Corollary 10.1 CBMP can be solved in 0{M{n,m) +nm) time, where n := \ V'\ and m := \E'\. 

When the input functions u, b in CBMP are small enough, the transit capacities of nodes in 
become small as well. Then the application of the shortest blocking IS-flow method may result 
in a competitive or even faster algorithm for CBMP. Let the capacities of all edges of G' be ones. 
We have A(A^) = 0{m) in general, and A(A) = 0{n) if b is all-unit. Then Corollarv 19.21 vields 
the same time bounds as in |121 125j for the corresponding cases. 

Corollary 10.2 SBFM (with the fast implementation of a phase as in Section\^ solves the max- 
imum degree- constrained subgraph (or b-f actor) problem in 0{m^/'^) time and solves the maximum 
matching problem in a general graph in 0{\/nm) time. 

Feder and Motwani ^ elaborated a clique compression technique and used it to improve the 
0{^/nm) bound for the maximum bipartite matching problem to 0(-^/nm log(n^/m)/ log n). We 
explain how to apply a similar approach to a special case of MSFP, lowering the bound for dense 
nonbipartite graphs. We need a brief review of the method in [^. 

Let H = {X,Y,E) be a bipartite digraph, where E (1 X x Y , \X\ = \Y\ = n and \E\ = m. 
A (bipartite) clique is a complete bipartite subgraph {A,B,A x B) of H, denoted by G{A,B). 
Define the size s(C) of C = G{A,B) to be |^| + \B\. A clique partition of is a collection C of 
cliques whose arc sets form a partition of E; the size s{C) of C is the sum of sizes of its members. 

Let a constant < S < 1/2 be fixed. Then a clique G{A,B) of H is called a 6-clique if 
1^1 = [n^"''] and \B\ = [51ogn/log(2n^/m)J. It is shown in P| that a 5-clique exists. 

The clique partition algorithm in |^ finds a <5-clique Ci in the initial graph Hi = (A, Y, E =: 
El ) and deletes the arcs of Gi , obtaining the next graph H2 = (A, y, £^2 ) • Then it finds a 5-clique 
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6*2 (concerning the number of arcs of H2) and delete the arcs of C2 from if2, and so on while the 
number of arcs of the current graph is at least 2n^~^ and the y-part of a 5-clique is nonempty. 
The remaining arcs are partitioned into cliques consisting of a single arc each. So the cliques 
Ci extracted during the algorithm form a clique partition. The running time of the algorithm is 
estimated as the sum of bounds T(Cj) on the time to extract the cliques Cj plus a time bound r' 
to maintain a certain data structure (so-called neighborhood trees). One shows that 

(29) the algorithm runs in 0{y/nm(3) time and finds a clique partition C oi H such that 
s(C) = 0{mP), where /3 := log(n^/m)/logn. 

Suppose we wish to find a maximum matching in a bipartite graph or, equivalently, to find a 
maximum integer flow from s to t in a digraph G with unit arc capacities, node set X UY U {s,t} 
and arc set E U (s x X) U (y x t), where E X x Y . One may assume \X\ = \Y\ = n. Using 
the above algorithm, form a clique partition C as in (|29() for {X,Y,E). Transform each clique 
C{A, B) in C into a star by replacing its arcs by a node z, arcs (x, z) for aW x ^ A and arcs (z, y) 
for all y E -B. There is a natural one-to-one correspondence between the s to t paths in G and 
those in the resulting graph G* , and the problem for G* is equivalent to that for G. Compared 
with G, the graph G* has \C\ additional nodes but the number m* of its arcs becomes 2n + s{C), 
or 0{m(3). Given a flow in G* , any (simple) augmenting path of length q meets exactly [q — l)/2 
nodes va. XiJY , and these nodes have unit transit capacities. This implies that Dinits' algorithm 
has 0{y/ri) phases (arguing as in [Hl[^). Since each phase takes 0{m*) time, the whole algorithm 
runs in 0{^/nmj3) time, as desired. 

Now suppose H = (X, Y, E, a) is a skew-symmetric bipartite graph without parallel arcs, 
where the sets X and Y are symmetric each other. We modify the above method as follows. 
Note that any two symmetric cliques in H are disjoint (otherwise some x G X is adjacent to cr(x), 
implying the existence of two arcs from x to o"(x)). We call a clique partition C symmetric if C G C 
implies cr(G) S C. An iteration of the symmetric clique partition algorithm works as in the previous 
algorithm, searching for a (^-clique G' in the current H, but then deletes the arcs of both G' and 
cr{G'). Let the algorithm construct a partition C consisting of cliques G[,a{G[), . . . ,G'^,a{Gl.) 
obtained in this order. 

To estimate the size of C and the running time, imagine we would apply the previous algorithm 
to our H (ignoring the fact that H is skew-symmetric). Let the resulting partition C be formed 
by cliques Gi, . . . ,Gq (in this order). Note that for a bipartite graph with n nodes and m arcs, 
both the number e(C) of arcs of a (^-clique G and its size s{G) are computed uniquely, and these 
are monotone functions in m, as well as the above-mentioned time bound t(C) (indicated in jH]). 
Moreover, one can check that m' < m and e(C") 7^ imply m' — 2e(C") < m — e(C), where G' is 
a 5-clique in a graph with n nodes and m' arcs. Using these, we can conclude that r < q and that 
for i = l,...r, s{Gi) < s{Gi) and t{G'i) < T{G^). Then s(C') < 2s(C), implying s{C) = 0{mp), 
by (|29|) . Also the time of the modified algorithm is 0{y/nmP) (by (|29|) and by the fact that the 
above bound r' remains the same) plus the time needed to treat the symmetric cliques (t(C-), 
which is 0{m). 
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Finally, the graph H* obtained from H by transforming the cliques C(, cj(C(), . . . , C',, (y{C[) 
into stars has a naturally induced skew-symmetry. By the above argument, H* has 0{mP) arcs, 
and computing H* takes 0{y/nmf3) time. Apply such a transformation to the input graph of 
MSFP arising from an instance of the maximum matching problem. Arguing as in the bipartite 
matching case above and as in the proof of Lemma 18.41 we conclude with the following. 

Theorem 10.3 A maximum matching in a general graph with n nodes and m edges can he found 
in 0{y/nm\og{n? /m) / log n) time. 
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